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Abstract 

We have used 106 pb _1 of data collected with the Collider Detector at Fermilab to 
search for new particles decaying to dijets. We exclude at the 95% confidence level 
models containing the following new particles: axigluons and flavor universal colorons 
with mass between 200 and 980 GeV/c 2 , excited quarks with mass between 80 and 
570 GeV/c 2 and between 580 and 760 GeV/c 2 , color octet technirhos with mass be- 
tween 260 and 480 GeV/c 2 , W bosons with mass between 300 and 420 GeV/c 2 , and 
Eq diquarks with mass between 290 and 420 GeV/c 2 . 

PACS numbers: 13.85.Rm, 12.38.Qk, 14.70.Pw, 14.80.-j 



In this paper we extend a previous search ]I|] for narrow resonances in the dijet 
mass spectrum in pp collisions at a center-of-mass energy y^ =1.8 TeV. The previous 
search used 19 pb _1 of data collected in 1992-93 from run 1A of the Tevatron. This 
search uses 106 pb _1 of data collected in 1992-95 from both run 1A and run IB, and 
significantly extends our sensitivity to new particles. 

As before, we perform both a general search for narrow resonances and a specific 
search for axigluons 0, excited quarks ||, color octet technirhos j|, W, Z' 0, and 
E 6 diquarks ||. In addition, the flavor universal coloron [[?), a hypothesized massive 
gluon which couples equally to all quarks, is considered together with axigluons. The 
cross section for the coloron is always greater than or equal to that of the axigluon, 



so our axigluon limits will apply to the coloron as well. In models of supersymmetry 
in which the gluino is lighter than 5 GeV/c 2 , there can be dijet resonances resulting 
from squark decay || ||. We do not consider this model, since data from both our 
previous search and from a preliminary version of the present search has already been 
used to exclude a range of squark masses in the light gluino scenario || U . 

A detailed description of the Collider Detector at Fermilab (CDF) can be found 



elsewhere [01]]. We use a coordinate system with the z axis along the proton beam, 
transverse coordinate perpendicular to the beam, azimuthal angle 0, polar angle 8, 
and pseudorapidity r\ = — lntan(#/2). Jets are reconstructed as localized energy de- 
positions in the CDF calorimeters that are arranged in a projective tower geometry. 
The jet energy E is defined as the scalar sum of the calorimeter tower energies inside 



a cone of radius R = J(Ar]) 2 + (A0) 2 = 0.7, centered on the jet direction. The jet 
momentum P is the corresponding vector sum: P = J2 EiUi with Ui being the unit 
vector pointing from the interaction point to the energy deposition E( inside the same 
cone. E and P are corrected for calorimeter non-linearities, energy lost in uninstru- 
mented regions of the detector and outside the clustering cone, and energy gained 
from the underlying event and multiple pp interactions. The jet energy corrections 
increase the jet energies on average by roughly 24% (19%) for 50 GeV (500 GeV) 
jets. Full details of jet reconstruction and jet energy corrections at CDF can be found 



elsewhere 11 



We define the dijet system as the two jets with the highest transverse momentum 



in an event (leading jets) and define the dijet mass m = y {E\ + E 2 ) 2 — (Pi + P2) 2 . 



The dijet mass resolution is approximately 10% for dijet mass above 150 GeV/c 2 . 
Our data sample was obtained using four triggers that required at least one jet with 
uncorrected cluster transverse energies of 20, 50, 70 and 100 GeV, respectively. After 
jet energy corrections these trigger samples were used to measure the dijet mass 
spectrum above 180, 241, 292 and 388 GeV/c 2 , respectively. At these mass thresholds 
the trigger efficiencies were greater than 95%. The four data samples corresponded to 
integrated luminosities of 0.126, 2.84, 14.1 and 106 pb _1 after prescaling. Offline we 
required that both jets have pseudorapidity \rj\ < 2 and a scattering angle in the dijet 
center-of-mass frame | cos 6**| = | tanh[(r/i — r] 2 )/2]| < 2/3. The cos 6* requirement 
provides uniform acceptance as a function of mass and reduces the QCD background 
which peaks at | COS6 1 *! = 1. To utilize the projective nature of the calorimeter towers, 
the z position of the event vertex was required to be within 60 cm of the center of the 
detector; this cut removed 7% of the events. Backgrounds from cosmic-rays, beam 
halo, and detector noise were removed with the cuts reported previously |jl|, and 
residual backgrounds were removed by requiring that the total observed energy be 
less than 2 TeV. 

In Fig. |I] we present the inclusive dijet mass distribution for pp — > 2 jets + X, 
where X can be anything including additional jets. The dijet mass distribution has 
been corrected for trigger and z vertex inefficiencies. We plot the differential cross 
section versus the mean dijet mass in bins of width approximately equal to the dijet 
mass resolution (RMS~ 10%). The data are compared to a QCD prediction from 
the PYTHIA Monte Carlo \Y1 and a simulation of the CDF detector. The cross 
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section predicted by the QCD simulation, using CTEQ2L parton distributions |T3 
and a renormalization scale fi = Pt, is normalized to the data in the first 6 bins 
(180 < m < 321 GeV/c 2 ) by dividing the simulation by a factor of 0.66. In Fig. [I] 
the horizontal lines on the data points indicate the bin width, the same width in data 
and simulation. The points are plotted at the mean mass, calculated independently 
for data and simulation. 

We note that the data is above the QCD simulation at high mass. In a previous 



paper Jl4|], we reported a similar effect in the fully corrected inclusive jet transverse 
energy distribution compared to an O(a^) parton level QCD calculation. Unlike 
the inclusive jet analysis, here we do not deconvolute the mass distribution for the 
effects of detector resolution, and instead compare the data directly to QCD plus a 
CDF detector simulation. In our previous dijet mass search the excess was not as 
noticeable because we normalized the simulation to the data on average, while here 



we normalize to the low mass end as described above. In another paper |15[ we have 
studied the dijet angular distributions and find them to be in good agreement with 
QCD in all regions, including at high mass. The source of the high dijet mass and 
high jet transverse energy excess is not yet fully understood. Candidate explanations 
within the Standard Model include a larger than expected gluon distribution of the 



proton [16] or large QCD corrections from resummation |17]]. As in our previous 
search |1[, we do not use QCD calculations to determine the background to new 
particles, but merely use the data itself to fit for the background. 

To search for resonances we fit the data with the parameterization da /dm = 
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A(l — m/ \fs + Cm 2 / ' s) jm F with parameters A, C, N and P. In the run 1A search [jl[ 
the term Cm 2 / s was not used because fewer parameters were needed to fit the lower 
statistics sample. With the higher statistics in this sample the extra term Cm 2 / s 
was needed to obtain an acceptable fit. This parameterization gives an adequate 
description of both the observed distribution (x 2 /DF = 1.49) and the QCD prediction 
(x 2 /DF= 0.85). Figure [l] shows the background fit on a logarithmic scale, and Fig. ||] 
shows the fractional difference between the data and background fit on a linear scale. 
Figures [L] and || also show the predicted line shape for excited quarks (q*) using 



the PYTHIA Monte Carlo |T2[ and a CDF detector simulation. If excited quarks were 
produced in pp collisions, their production and decay to dijets would proceed via the 
process qg — > q* — » qg. The mass resolution is dominated by a Gaussian distribution 
(RMS~ 10%) from jet energy resolution and a long tail towards low mass from QCD 
radiation. Since the natural width of a q* is significantly smaller than the measured 
width, the q* mass resonance curves in Figs. |l| and |] were used to model the shape 
of all narrow resonances decaying to dijets. 

There is no statistically significant evidence for a dijet mass resonance, which 
should appear in at least two neighboring bins above the background fit. We note 
that in the region of 550 GeV/c 2 there is a single bin which is 2.6 standard deviations 
above the fit; however, this region is not well fit by a new resonance because the 
number of events in neighboring bins is too low. When we fit the data to both a 550 
GeV/c 2 resonance and a smooth background we find that the upward fluctuation in 
the data is significantly narrower than expected for a resonance. 



Systematic uncertainties on the cross section for observing a new particle in the 
CDF detector are shown in Fig. [2| Each systematic uncertainty on the fitted sig- 
nal cross section was determined by varying the source of uncertainty by ±ler and 
refitting. In decreasing order of importance the sources of uncertainty are the 5% 
jet energy scale uncertainty, low mass data, the background parameterization, QCD 
radiation's effect on the mass resonance line shape, trigger efficiency, jet energy reso- 
lution, relative jet energy corrections between different parts of the CDF calorimeter, 
energy scale of run 1A with respect to run IB, luminosity and efficiency. For example, 
at 600 GeV/c 2 reducing the jet energy by 5% centers the resonance on an upward 
fluctuation, and increases the fitted signal by 225%. The low mass data uncertainty, 
listed above, is because the background fit gets significantly worse when data between 
150 and 180 GeV/c 2 are included. The larger number of interactions per crossing in 
run IB increases the uncertainty on the lower mass data, so we start the mass distri- 
bution at 180 GeV/c 2 . However, since this mass range was included in run 1A, the 
effect of adding the low mass data is included as a systematic for run 1A plus run IB. 

The total systematic uncertainty was found by adding the individual sources in 
quadrature. In this analysis the relative systematic error is larger than it was in the 
previous analysis: the total run 1A and IB systematics range from 40% to 300% of 
the cross section while the run 1A systematics ranged from 30% to 120%. This is not 
because the absolute systematics have significantly increased, but instead because the 
size of the signal we are statistically sensitive to has decreased by over a factor of 
two, so now the systematics have a larger relative effect. This is particularly true at 



masses near upward fluctuations in the data. 

In the absence of conclusive evidence for new physics we proceed to set upper 
limits on the cross section for new particles. For each value of new particle mass in 
50 GeV/c 2 steps from 200 to 1150 GeV/c 2 , we perform a binned maximum likelihood 
fit of the data to the background parameterization and the mass resonance shape. 
We convolute each of the 20 likelihood distributions with the corresponding total 
Gaussian systematic uncertainty, and find the 95% confidence level (CL) upper limit 
presented in Table I. 

In Fig. [| we plot our measured upper limit on the cross section times branching 
ratio for a new particle decaying to dijets as a function of new particle mass in 50 
GeV/c 2 steps. The points are connected by a smooth curve, which is an estimate of 
the upper limit in between the measured points. The limit is compared to lowest order 
theoretical predictions for the cross section times branching ratio for new particles 
decaying to dijets IJ. New particle decay angular distributions are included in the 
calculations, and we required \rj\ < 2 and |cos#*| < 2/3 for all predictions. For 
axigluons (or flavor universal colorons) we exclude the mass range 200 < Ma < 980 
GeV/c 2 , extending the previous CDF exclusions of 120 < M A < 870 GeV/c 2 0. 
For excited quarks we exclude the mass ranges 200 < M* < 520 and 580 < M* < 
760 GeV/c 2 , significantly extending the previous CDF exclusion of 80 < M* < 570 
GeV/c 2 [[[], 0. The DO collaboration has performed a preliminary search for excited 
quarks and exclude the mass range 200 < M* < 720 GeV/c 2 [III. These exclusions 



are for Standard Model couplings (/ = /' = f s — !)• For smaller couplings, the new 
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excluded region in the coupling [[J vs. mass plane is shown in Fig. [| compared to 
previous excluded regions. For color octet technirhos (pt) we exclude the mass range 
260 < M PT < 470 GeV/c 2 , extending to lower mass the previous CDF exclusion of 
320 < M PT < 480 GeV/c 2 |l|. For the first time we exclude the hadronic decays of 
the new gauge boson W in the mass range 300 < Mw r < 420 GeV/c 2 . Also for the 
first time we exclude E$ diquarks in the mass range 290 < Me s < 420 GeV/c 2 . The 
cross section for hadronic decays of Z' is too small to exclude. 

In conclusion, the measured dijet mass spectrum does not contain evidence for 
a mass peak from a new particle resonance. We have presented model independent 
limits on the cross section for a narrow resonance, and set specific mass limits on ax- 
igluons, flavor universal colorons, excited quarks, color octet technirhos, new charged 
gauge bosons, and E§ diquarks. 
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Table I: As a function of new particle mass we list our 95% CL upper limit on 
cross section times branching ratio for narrow resonances decaying to dijets. The 
limit applies to the kinematic range where both jets have pseudorapidity \q\ < 2.0 
and where the dijet system satisfies | cos#*| < 2/3. 
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Figure 1: The dijet mass distribution (circles) compared to a QCD simulation (boxes) 
and fit to a smooth parameterization (solid curve). Also shown are simulations of ex- 
cited quark signals in the CDF detector (dashed curves). In the data and simulations 
we require that both jets have pseudorapidity \v\ < 2.0 and that the dijet system 
satisfies | cos 0*1 < 2/3. 
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Figure 2: The fractional difference between the dijet mass distribution (points) and a 
smooth background fit (solid line) is compared to simulations of excited quark signals 
in the CDF detector (dashed curves). The inset shows the systematic uncertainty for 
a new particle signal (see text). 
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Figure 3: The upper limit on the cross section times branching ratio for new particles 

decaying to dijets (points) is compared to theoretical predictions for axigluons [0, 

flavor universal colorons J7|, excited quarks ||, color octet technirhos ||, new gauge 

bosons W and Z' ||, and Eq diquarks 0. The limit and theory curves require that 

both jets have pseudorapidity \r}\ < 2.0 and that the dijet system satisfies | cos^*! < 

2/3. 

17 



CO 



CO 



0.9 



Q. 

O 0.8 

o 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



LEP 




100 



200 



300 



400 



500 



600 



700 



800 



Excited Quark Mass (GeV/c ) 



Figure 4: The region of the coupling vs. mass plane excluded by previous CDF 
measurements |L[ |18|] in the q* — > 97 and g* — ► qW channels (clear hatched region) 
and q* —>■ qg channels (shaded hatched region on left) in run 1A is extended by this 
q* — > qg search in run 1A plus run IB (shaded hatched region on right). The CDF 
excluded regions are compared to the regions excluded by LEP I and UA2 (shaded 
regions) [gjfj. 1S 



